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ABSTRACT 

In the nonlinear regime of cosmic clustering, the mass density field of the 
cosmic baryon fluid is highly non-Gaussian. It shows different dynamical behav- 
ior from collisionless dark matter. Nevertheless, the evolved field of baryon fluid 
is scale-covariant in the range from the Jeans length to a few ten h _1 Mpc, in 
which the dynamical equations and initial perturbations are scale free. We show 
that in the scale-free range, the non-Gaussian features of the cosmic baryon fluid, 
governed by the Navier-Stokes equation in an expanding universe, can be well 
described by a log-Poisson hierarchical cascade. The log-Poisson scheme is a 
random multiplicative process (RMP), which causes non-Gaussianity and inter- 
mittency even when the original field is Gaussian. The log-Poisson RMP contains 
two dimensionless parameters: (3 for the intermittency and 7 for the most singu- 
lar structure. All the predictions given by the log-Poisson RMP model, including 
the hierarchical relation, the order dependence of the intermittent exponent, the 
moments, and the scale-scale correlation, are in good agreement with the results 
given by hydrodynamic simulations of the standard cold dark matter model. The 
intermittent parameter f3 decreases slightly at low redshift and indicates that the 
density field of baryon fluid contains more singular structures at lower redshifts. 
The applicability of the model is addressed. 

Subject headings: cosmology: theory - large-scale structure of universe 



1. Introduction. 

In the universe, about 72% of the energy density is in the form of dark energy, 24% cold 
dark matter, and a small fraction, 4% baryon matter. Since the dark energy is assumed to be 
spatially uniform, the dynamics of the clustering of cosmic baryon fluid should be dominated 
by the underlying gravitational potential of dark matter. However, it has already been 
recognized in the early study of cosmic structure formation that in the nonlinear regime the 
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dynamical behavior of the cosmic baryon fluid, or of the intergalactic medium (IGM), doesn't 
always follow the collisionless dark matter. Although the cosmic baryon fluid is passive 
substance in comparing with dark matter, it statistically decouples from the underlying dark 
matter field in the non-linear evolutionary stage. In the scale free range, the cosmic baryon 
fluid, as a Navier-Stokes fluid in the expanding universe, is similar to the fluid being moved 
by inertia, and should show some features as the turbulence in inertial range (Shandarin and 
Zeldovich 1989). 

Later, it was found that the dynamical equations of the velocity fields of cosmic matter 
essentially are a variant of the random-force-driven Burgers' equation (Gurbatov et al. 1989; 
Berera & Fang 1994). For baryon fluid, it is a Burgers' equation driven by the random force 
of the gravity of dark matter (Jones 1999; Matarrese & Mohayaee 2002). Burgers' fluid will 
show highly non-Gaussian features due to the development of Bergers' turbulence when the 
Reynolds number is large enough (Polyakov 1995; Lassig 2000; Bee & Frisch 2000; Davoudi 
et al. 2001). In this state, the Burgers' fluid consists of shock waves in low as well as in high 
density regions, and therefore, non-Gaussianity can be seen in low as well as in high density 
regions. 

This property has received supports from the absorption spectra of QSOs, which is 
caused by the IGM with moderate mass density. For example, the Lya transmitted flux in 
the absorption spectra of QSOs is found to be significantly intermittent and its probability 
distribution functions (PDF) are remarkably long tailed (Jamkhedkar et al. 2000; Pando et 
al. 2002; Feng et al. 2003). The H I and He II Lya absorption lines of QSO HE2347 can 
not be explained by thermal broadening, but consistent with turbulence broadening (Zheng 
et al. 2004; Liu et al. 2006). Moreover, samples produced by cosmological hydrodynamic 
simulations also reveal the statistical decoupling of baryon matter from dark matter and the 
non-Gaussian features of Bergers' turbulence (He et al. 2004, 2005; Kim et al. 2005). 

One of the latest developments in this direction is that the random velocity fields of the 
cosmic baryon fluid are found to be extremely well described by She-Leveque's (SL) scaling 
formula (She & Leveque 1994) in the scale range from the Jeans length to larger than 10 h _1 
Mpc (He et al. 2006). The SL scaling formula is believed to characterize the scaling hierarchy 
of the non-linear evolution of Navior-Stokes fluid, like fully developed turbulence (e.g., Frish 
1995). Therefore, the dynamical features of the cosmic baryon fluid in the nonlinear regime 
are similar to the fully developed turbulence: it is of scaling hierarchy. It is interesting to 
note that the SL formula has also been successfully applied to describe the mass fields of gas 
on interstellar scales (Boldyrev et al. 2002; Padoan et al. 2003). 

The SL formula actually is originated from a cascade of log-Poisson random multiplica- 
tive process (RMP), which is related to the hidden symmetry of the Navier-Stokes equations 
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(Dubrulle 1994; She & Waymire 1995; Benzi et al. 1996). This motivate us to investigate 
whether the clustering behavior of the mass density field of the cosmic baryon fluid can 
be described by the log-Poisson RMP. Theoretically, this is not trivial, because the cosmic 
baryon fluid is compressible and dominated by the gravitational field of dark matter, while 
the SL scaling formula originally was proposed to describe the velocity field of incompressible 
fluid. 

The paper is organized as follows. §2 describes the model of the log-Poisson RMP for 
the nonlinear evolution of the cosmic baryon fluid. The predictions and its tests of the log- 
Poisson RMP model are presented in §3. §4 presents briefly the redshift-evolution of the 
coefficients of the log-Poisson model. Discussion and conclusion are given in §5. 



2. Log-Poisson RMP model 

2.1. The log-Poisson hierarchy 

The clustering and non-Gaussianity of the cosmic mass density and velocity fields are 
usually measured by two and multiple point correlation functions. To reveal the features 
of the scaling hierarchy of the mass density field, however, it is more effective to use the 
structure function defined by 

S p (r) = (1**1*), (1) 

where 5p r = p(x + r) — p(x), r = |r|, p is the order of statistics, and the average (...) is taken 
over the ensemble of density field. For statistically isotropic and homogenous random field 
p(x), S p (r) depends only on r. 

The difference between the correlation function and structure function has been analyzed 
in detail by Monin & Yaglom (1975). The variable 5p r = p(x + r) — p(x) is not 5p(x) = 
p(x) — p, p being the mean of density; the variable 5p(x) can be larger than p, but cannot be 
less than — p, and therefore, for a nonlinear field, the distribution of 5p(x) generally is skew; 
however, the distribution of Sp r is symmetric with respect to positive and negative 5p r if the 
field is statistically uniform. 

In the scale-free range of the fluid, the structure function as a function of r can be 
expressed as a power law 

S p (r) oc r ?(p) . (2) 

For fully developed turbulence of Navier-Stokes fluid, £(p) is a nonlinear function of p, i.e. 
the mass field is intermittent, and £(p) is called intermittent exponent (Frisch 1995). Since 
the pioneer work of Kolmogorov (1941), it is believed that the relation of £(p) vs. p is 
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related to the scale-covariance of the dynamical equations and initial conditions. Since then 
many hierarchy models for interpreting have been proposed. Finally the best model 
is given by the SL scaling formula (She & Leveque 1994). It has been shown that the SL 
formula is yielded from the Log-Poisson hierarchy process, which is related to the so-called 
generalized scale covariance of the Navier-Stokes equations (Dubrulle 1994). Therefore, one 
may expect that the statistical behavior of the mass field of cosmic baryon matter would 
also be interpreted by the log-Poisson random multiplicative processes (RMP) . 

The log-Poisson RMP assumes that, in the scale-free range, the variables \Sp r \ on dif- 
ferent scales r are related from each other by a statistically hierarchy relation given by 

\5p r2 \ = W rir2 \5p ri \, (3) 

where 

W rir2 = /r(n/r 2 r, (4) 

which describes how the fluctuation \5p ri \ on the larger scale r\ related to fluctuations \Sp r2 \ 
on the smaller scale r 2 . In eq.(4), m is a Poisson random variable with the PDF 

P(m) = exp(-A rir2 )A™ r2 /m!. (5) 

To insure the normalization (W rir2 ) = 1, where (...) is over m, the mean A rir2 of the Poisson 
distribution should be 

A rir2 = 7 [ln(r 1 /r 2 )]/(l-/3). (6) 

It is enough to consider only \5p r \, as the distribution of positive and negative 5p r is sym- 
metric. 

The log-Poisson model of equation (3) depends only on the ratio r\jri, thus, it is scale 
invariant. The model is determined by two dimensionless positive parameters: (3 and 7, of 
which the physical meaning will be given below. Equation (3) relates Sp r on different scales 
by multiplying a random factor W, and therefore, it is a random multiplicative process 
(RMP), which generally yields a non-Gaussian field even when the field originally to be 
Gaussian (Pando et al. 1998). For a Gaussian field, variables Sp ri and 5p T2 are statistically 
independent and it requires (5 — > 1 in equation (4). 



2.2. Parameter 7 and singular structures 

With the log-Poisson model, the intermittent exponent is given by (see Appendix) 

£(p) = -7[p-(l-/3 p )/(l-/3)]- (7) 
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For a Gaussian field with scale-free power spectrum, we have oc p. This indicates again 
that a Gaussian field requires (3 — > 1. 

Considering a density field containing singular structures at positions x i; we have p(x) oc 
|x — Xi|~ Q near x ; with a > 0. The variables \Sp r \ p near x ; should be ~ |r|~ ap , and therefore 
\Sp r \ p+1 /\Sp r \ p ~ \r\~ a . Thus, to pick up the singular structures, we define a statistical tool 

as 

F p (r) = S p+1 (r)/S p (r). (8) 

For higher p, singular structures have larger contributions to S p (r), while for lower p, weak- 
clustering structures have larger contributions to S p (r). Therefore F p (r) measures clustering 
structures, which are dominant for the p-order statistics. Obviously, when p — > oo, F p {r) 
should be dominated by the singular structures, i.e., oc r~ a . 

On the other hand, from equations (2), (7), and (8), one finds 

F p (r) oc r'^ 1 -^. (9) 

Since (3 < 1, we have 

F = lim ^ 6pr ^ oc r" 7 (10) 

(\s Pr \p) ^ r ■ [W) 

Therefore, the parameter 7 of the log-Poisson RMP is actually the power-law index of the 
mass profile of the most singular structures. It should be pointed out that the word "singular" 
is applicable only asymptotically, because we cannot let r — > to pick up the singular 
structure, as r should be in the scale-free range. 



2.3. Parameter (3 and intermittency 

As mentioned in last subsection, when /3 — > 1, the field would be Gaussian, and therefore, 
equation (9) implies that for a Gaussian field, F p (r) would be r-independent. On the other 
hand, for a utmost intermittent field, which consists only the most singular structures, and 
the fluctuations between the singular structures are zero, F p (r) should be equal to r~ 7 
regardless of p. From equation (9), the utmost intermittent field should have (3 = 0. Thus, 
the parameter (3 is to measure the level of intermittency: non-intermittency corresponds to 
(3 = 1, and the strongest intermittency corresponds to (3 = 0. 

The meaning of (3 can also be seen with the following hierarchical relation of F p (r) 

Foo(r) 



F P+ i(r) 
Foo(r) 



VP 
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which can be derived from equations (9) and (10). Equation (11) is invariant with respect 
to a translation in p. In deriving eq.(ll), we assume that the proportional coefficient of 
equation (9) is p-independent. We will show that this assumption is correct. 

As mentioned in §2.2, the F p (r) measures the clustered structures dominating the p 
order statistics. The smaller the p, the larger the contribution of weak-clustering structures 
to F p (r). Therefore, equation (10) describes the hierarchical relation between the stronger (or 
high p) and weaker (or low p) clustering. In the scale-free range where F p+1 {r) / F^ir) < 1, 
we have F p (r)/F 00 (r) < F p+1 (r)/F 00 (r) if /3 < 1. That is, for an intermittent field, weak 
clustering structures are strongly suppressed with respect to the most singular structures; 
the smaller the (3, the stronger the suppression of weak clustering structures. 



The samples for testing the model of §2 are similar to that used in He et al. (2006), which 
are given by a hybrid hydrodynamic/iV-body simulation, consisting of the WENO algorithm 
for baryon fluid and iV-body simulation for particles of dark matter (Feng et al. 2004). We 
now produce samples of 50 h~ l Mpc box, 768 3 grid, and the cosmo logical parameters are 
taken from the results of WMAP (Spergel et al. 2006). The samples output at redshifts 0, 
1, 2, 3, and 4. We randomly sample 10,000 one-dimensional sub-samples at each redshift. 

In order to have a complete description of the density field and avoid false statistical 
correlation, the variable 5p r should be given by a proper decomposition of the field p x . 
We will use the decomposition of discrete wavelet transform (DWT), which is found to be 
effective to describe turbulence (e.g., Farge 1992). With the DWT, the variables of mass 
density field is given by 



where ^j,i(x) is the base of discrete wavelet transform (e.g., Fang & Thews 1998). For a 
one-dimensional sample of length L, the scale index j is related to the scale r by r = L/2 J 
and the position index / is for the cell at x = lL/2 j to (/ + l)L/2 j . We will use the Harr 
wavelet to do the calculation below. We also repeat the calculations with wavelet Daubechies 
4. The non-Gaussian statistical features given by Daubechies 4 are the same as that of Haar 
wavelet. 



3. Non-Gaussianity of the cosmic baryon fluid 




(12) 



-7- 



10' 




8.3 h~'Mpc 

8. 3, fitted Gauss. 

4.2 

2.1 




10 



5 





nl5p r 



5 



10 



Fig. 1. — PDFs of the density difference ln|<5p r | with r = 2.1, 4.2 and 8.3 h 1 Mpc for 
samples at z — 0. A fitted Gaussian PDF of In \ 5p r \ is also shown. 



We first show the basic statistical deviation of the baryon mass density field from a 
Gaussian field. Figure 1 gives the PDF of the density difference variables, p(ln \5p r \)d\n \8p r \, 
for the cosmic baryon fluid sample at z = on scales r = 2.1, 4.3, 8.3 hr l Mpc. A fitted 
Gaussian to the PDF on scale 8.3 ft, _1 Mpc is also shown as the thick solid line in Figure 1. 
Since the Gaussian fitting is for In \Sp r \, the fitted curve actually is a lognormal distribution 
for \5p r \. It shows clearly that on all the scales, the PDFs of \Sp r \ are non-Gaussian and 
have a longer tail than the lognormal distribution. 

The long tailed events can be effectively described by the p-dependence of the structure 
function S p (r). Figure 2 shows S p (r) as a function of r for p = 0.5 to 4 for the cosmic baryon 
fluid sample at z — 0. For all p the structure function \nS p (r) of fig. 2 can be well fitted 
by a straight line in the scale range of 2 < r < 16 h _1 Mpc. The data points at 1 h _1 
Mpc are slightly deviating from the straight line given by the fitting over 2 < r < 16 h _1 
Mpc, because 1 h" 1 Mpc is already close to the Jeans length. We only focus on the range of 
2 < r < 16 h _1 Mpc below. The upper limit 16 h _1 Mpc actually is from the finite size of 
the simulation box. 

For p > 1, the structure function decreases when the scale r increases from 2 to 16 h -1 
Mpc, while for p < 1 it increases with the increase of scale r. This requires < for 
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Fig. 2. — Structure functions S p (r) vs. r in the range 1 < r < 16 h" 1 Mpc. p is equal to 
0.5 x n and n — 1, 2. ..8 from bottom to top. 

p > 1, and > for p < 1. Therefore, the intermittent exponent can't be fitted by 
P- This is once again to show that the field is highly non-Gaussian. On the other 
hand, equation (7) does show that can have different signs for p > 1 and p < 1. 



3.2. Hierarchical relation and parameter (3 

We now test the hierarchical relation (11), which can be rewritten as 

lnF p+1 (r) =p\nF p (r) + A(r), (13) 

where A(r) = (1 — (3) lnF^r) depends only on r. Quantity A(r) may also depend on p if 
the proportional coefficient of relation (9) is p-dependent. 

Equation (13) requires that lnF p+1 (r) vs. \nF p (r) should be a straight line for a given 
scale r, and the slope (3 should be the same for all the scales r. Figure 3 presents the relations 
of lnF p+ i(r) vs. In F p (r) of the samples at z = for the scales r = 2.1, 4.2, 8.3 and 16.7 
h _1 Mpc and for p — 1, 1.5, 2, and 2.5. It shows that all the relations of lnF p+1 vs. lnF p at 
different r can be well fitted by straight lines with slope (3 = 0.28 ± 0.02. 

Figure 3 also shows that A(r) depends only on r, but not on p. This is consistent with 
the assumption of the p- independence of the proportional coefficient of relation (9). If A(r) 
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Fig. 3. — lnF p+ i vs. lnF p . r is equal to 2.1, 4.2, 8.3, and 16.7 h" 1 Mpc for lines from top to 
bottom. In each lines, the four data points correspond to p — 1, 1.5, 2, and 2.5 . 

is p-independent, we can find the following relation from eq.(13) 

HF p+1 (r)/F 3 (r)} = 0]n[F p (r)/F 2 (r)]. (14) 

This is, for all r and p, ln[F p+ i(r)/F 3 (r)] vs. ln[F p (r) / F 2 (r)] should be on a straight line. 
The relation of equation (14) is tested in Figure 4. All data points can indeed be fitted by 
a straight line with slope 0.28. This is the hierarchy of the density fluctuations between 
different order p on various scale r. 

3.3. Intermittent exponent and parameter 7 

The intermittent exponent £(p) as a function of p can be measured by fitting S p (r) 
(Figure 2) with a straight line of hi S p (r) = £(p) lnr + const for each p. The measured £(p) 
for the sample at z = are shown in Figure 5. The error bars are the variance of £(p) over 
100 samples, each of which contains 100 one-dimensional sub-samples. 

Equation (7) shows that the shape of as a function of p depends only on parameter 
f3, while parameter 7 gives the overall amplitude of the curve Since (5 is already 

determined in the last section, we can determine the parameter 7 by fitting Equation (7) to 
the amplitude of the measured The best fitting result is 7 = 0.91. The fitted curve 

are also shown in Figure 5. It shows that the feature of the intermittent exponent 
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Fig. A.— ln[F p+ i(r)/F 3 (r)] vs. ln[F p (r) / F 2 (r)] for data points of r =2.1, 4.2, 8.3 and 16.7 
h" 1 Mpc and p = 1, 1.5, 2, and 2.5. 

of the cosmic baryon fluid at z = in the range of 0.5 < p < 6 can be well reproduced with 
the log-Poisson model with parameters (3 = 0.28 and 7 = 0.91. 

One can make a further test of parameter 7 from equation (9). Since (3 = 0.28, we have 
(0.28) 3 ~ 0.02 < 1, and then, equation (9) yields 



Equation (15) requires that the relations of In F p (r) vs. lnr should be straight lines for all 
orders p > 3 with the same slope of 7. Figure 6 presents the relation between lnF p (r) and 
lnr, which can be fitted by straight lines in the scale range of 2 < r < 16 h _1 Mpc. The slope 
of the lines with p > 3 are 0.88 ± 0.06, consistent with the value of 7 = 0.91 determined 
from (Figure 5). 



With the determined parameters (3 and 7, we can predict statistical properties of the 
cosmic baryon fluid without other free parameters. As the first one, we consider the ratio 
between the high order and 2nd order moments, {5 p^ p ) / (5 p 2 r .) p . When p — 2, it is kurtosis, 
which is a popular tool to detect non-Gaussianity. For a Gaussian field, the ratio should be 
a constant, independent of r. The 2nd moment (5p 2 ) actually is the power spectrum of the 



In F p (r) 



~ —7 lnr + const, if p > 3. 



(15) 
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Fig. 5. — Intermittent exponent The solid line is given by equation (7) with (3 = 0.28 

and 7 = 0.91. The data points are from the fitting to the structure function S p (r) with 
equation (7). The error bars are the variance of over 100 samples, each of which 
contains of 100 one-dimensional sub-samples. 



mass density field (Fang & Feng 2000). For the log-Poisson model we have (see Appendix) 

In 7/5/ = K P In r + const (16) 
{5p 2 r )P 

with 

Kp = (17) 

That is, ln((5p 2p ) / (5p 2 .) p ) is linearly dependent on lnr (scale free) with the coefficient K p 
determined by (5 and 7. As expected, for Gaussian field (f3 — > 1), K p = 0, i.e., the ratio of 
moments is independent on lnr. 

Figure 7 shows the relation of ln((5p 2p ) / (5p 2 .) p ) vs. lnr for the sample at redshift z — 0. 
For clarity, we show only the results of p = 2 and 3, which correspond to the statistical order 
4 and 6. The errors are calculated as the variance over 100 samples, each of which contains 
100 lines. Since the error bars actually are very small in logarithm scale, one can not show 
them in Figure 7. The solid lines of Figure 7 are given by a least square fitting and have 
slopes 1.06 ± 0.06 and 2.10 ± 0.12, which are in agreement with the values 1.07 and 2.23 
calculated from equation (17) with (3 = 0.28 and 7 = 0.91. 
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Fig. 6. — lnFp(r) vs. lnr for p = 0.5 to 4 from bottom to top. 

3.5. Scale-scale correlation 

A powerful non-Gaussian detector is the so-called scale-scale correlation, which is defined 

as 

Obviously, for a Gaussian field, Cff r =1. It has been shown that one can construct a 
non-Gaussian field, which has identical first and second order statistics as a Gaussian field, 
but has strong scale-scale correlation (Pando et al 1998). 

It is especially important to measure the scale-scale correlation of cosmic baryon matter. 
The clustering of cosmic large scale structure in the nonlinear regime essentially is due to 
the interaction between Fourier modes on different scales (e.g., Peebles 1980). Therefore, 
cosmic clustering will definitely yield the transfer of the density perturbation powers between 
different scales, which leads to scale-scale correlation. Scale-scale correlation is also effective 
to distinguish various hierarchy cascade models (Pando et al. 1998). For instance, a Gaussian 
hierarchy cascade, like the model of Cole and Kaiser (1987), still yields = 1, while the 

so-called p-model and a-model yield C%f depending on both p and r (Greiner et al. 1996). 

If the ratio r 2 /Vi is fixed, the log-Poisson model predicts the scale-scale correlation to 
be (see Appendix) 

C™ 2 =B(r 2 / ri )4 i2p) ~ 2aP \ (19) 
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Fig. 7. — The ratio between high order and 2nd order moments (Sp 2p ) / (Spf,) p as a function 
of r for p = 2 and 3. The solid lines are given by the least square fitting, which yields slopes 
consistent with the ones calculated from equation (17) with (3 = 0.28 and 7 = 0.91. 

where the coefficient B(r 2 /r 1 ) depends only on the ratio r 2 /ri, as the log-Poisson model is 
invariant of the dilation. 

Thus, assuming rijr\ remains constant, the relationship of lnC^ p r vs. lnri should be 
a straight line with the slope £(2p) — 2£(p) = —7(1 — /3 P ) 2 /(1 — (3). The result is shown in 
Figure 8, in which we take ri/r 2 = 2 and p — 2, 3. When the number p = 3, the statistics of 
equation (19) actually is of 2p = 6 order. Figure 8 shows that for r > 4/i~ 1 Mpc, the scale- 
scale correlations can be well fitted by straight lines with slopes 1.11 ± 0.02 and 1.20 ± 0.01 
for p = 2 and 3, which are consistent with the values 1.10 and 1.20 calculated from eq.(19) 
with (3 = 0.28 and 7 = 0.91. The lower limit 4 hr l Mpc for the scale-scale correlation is 
higher than the lower limit 2 h _1 Mpc of the statistics in previous sections. It is because 
we take r 2 jr\ = 2, and the scale-scale correlation of 4 h^ 1 Mpc actually is the correlation 
between modes of 4 and 2 h _1 Mpc. 



4. Evolution of (3 and 7 

We repeat the similar analysis for samples at redshifts z — 1, 2, 3, and 4. The non- 
Gaussian features of all these samples can also be well explained with the log-Poisson model. 
The parameters (3 and 7 are listed in Table 1, which shows that both the parameters (3 and 7 
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Fig. 8. — Scale-scale correlation of the sample at z = for p = 2 and 3 with rijr\ = 2. 
The error bars are the variance over 100 samples, each of which contains of 100 lines. 

are increasing with redshift. The increase of /3 with redshift indicates that the intermittency 
is stronger at lower redshifts, and the fields at higher redshifts contain less singular structures 
than that at lower redshifts. 

On the other hand, the increase of 7 with redshift indicates that the singular feature is 
even stronger at higher redshift. This probably is because the baryon fluid is significantly 
heated by the Burgers' shocks at lower redshift (He et al. 2004) and leads to weaker singular 
structures. 



Table 1: (3 and 7 at different redshift z 



z 





1 


2 


3 


4 




0.28 


0.34 


0.34 


0.38 


0.43 


7 


0.91 


0.91 


1.0 


1.06 


1.16 
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5. Discussion and conclusion 

In the nonlinear regime of cosmic clustering, the dynamical behaviors of either dark mat- 
ter or baryon fluid are complicated. Nevertheless, it is generally believed that the evolution 
should be scale-covariant in the range where the dynamical equations and initial perturba- 
tions are scale-free. Therefore, hierarchical and universal scaling relations have been widely 
used to describe the nonlinear clustering. For instance, the hierarchical relations of irre- 
ducible correlation functions and the universal density profile of halos are successful in the 
description of the statistical features of massive halos of dark matter. 

However, it has already been recognized in the early study of cosmic structure formation 
that in the nonlinear regime the dynamical behavior of cosmic baryon doesn't always follow 
the collisionless dark matter. The non-Gaussianity of the mass and velocity fields of baryon 
fluid cannot be given by a similar mapping of the mass and velocity fields of dark matter. 
For instance, the halo model assumes that all mass fields are given by a superposition of the 
halos on various scales, and all non-Gaussian behaviors of the density field are described by 
the universal density profile (e.g., Cooray & Sheth 2002), this makes it difficult to explain 
the intermittency and the scale-scale correlation of the transmitted flux in the absorption 
spectra of QSOs. 

We show that the evolution of the cosmic baryon fluid, governed by the Navier-Stokes 
equation in an expanding universe, is also hierarchical in the scale range in which the dy- 
namical equations and initial perturbations are scale-free. The non-Gaussian behavior of 
the mass density field of baryon fluid can be well explained by the log-Poisson hierarchical 
cascade model. The SL formula and/or log-Poisson model are universal for the fully de- 
veloped turbulence of Navier-Stokes fluid in the scale-free range. Therefore, the result of 
this paper implies that, in the scale-free range, the cosmic baryon fluid reaches a statisti- 
cally quasi-steady state. For a fully developed turbulence, energy passes from large to the 
smallest eddies and finally dissipates into thermal motion, while the cosmic baryon fluid 
undergoes the evolution of clustering and finally falls into massive halos of dark matter to 
form structures, including light-emitting objects. Therefore, the log-Poisson model works on 
the scale range from the onset scale of the nonlinear evolution (a few tens of h _1 Mpc) to 
the dissipation scale, i.e., the Jeans length. 

In view of this picture, one can say that in the nonlinear regime, the statistical proper- 
ties of the cosmic baryon fluid are actually less dependent on the details of the dissipative 
processes. This property has already been noted in describing baryon matter by the Burg- 
ers' equation. Although the Burgers' equation contains a dissipative term, which leads to 
the formation of shocks and condense into luminous objects (Jones 1999), the self-similar 
properties of Burgers' turbulence actually depend very weakly on the dissipative term. 
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We now address the possible applications of the log-Poisson model. First, it is interesting 
to compare the log-Poisson model with the lognormal model, which assumes that the PDFs 
of the cosmic baryon matter is log-normal and no details of dissipative processes are needed 
(Bi & Davidsen 1997). The lognormal model is successful to explain some statistical features 
of Lyct forests and also predicts that the transmitted flux in the spectra of QSOs is non- 
Gaussian and intermittent. This result is qualitatively consistent with the observed data; 
however, non-Gaussian features given by the lognormal model do not quantitatively fit the 
observed data. For instance, the high order moment (§3.4) given by the lognormal model 
has a K p oc p(p — 1), while the data show K p oc — p 0A (p — 1) (Pando et al. 2002). The later 
is actually close to the log-Poisson model. Therefore, the higher order statistics of the Lya 
transmitted flux would be able to discriminate between the log-Poisson and the lognormal 
model. 

Second, recent studies have shown that the turbulence behavior of baryon gas can be 
detected by the Doppler-broadened spectral lines (Sunyaev et al. 2003; Lazarian & Pogosyan 
2006). Although these works focus on the turbulence of baryon gas in clusters, the result 
is still applicable, at least, for the warm-hot intergalactic medium (WHIM), which is shown 
to follow the evolution of Burgers' fluid on large scales (He et al. 2004, 2005). The last but 
not least, the polarization of CMB is dependent on the density of electrons, and therefore, 
the map of CMB polarization would provide a direct test on the non-Gaussian features of 
ionized gas when the data on small scales becomes available. 
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A. Log-Poisson model and intermittancy exponent 

In this appendix we give the details of deriving the statistical properties of the log- 
Poisson cascade model. Let us consider the log-Poisson model 

\Sp ri \ = W rori \5p ro \, (Al) 

where 

W rori = P m (r / ri )\ (A2) 
and m is a Poisson variables with probablity distribution function 

P(m)=exp(-A rori )A™ ri /m!. (A3) 



and 
Thus 



A rori = 7 HWn)]/(i-^). 

WU> = E(/ 3m ( r oM) 7 ) P exp(-A rori )A™ ri /m 

m 

= exp(-A ron )(r /r 1 )^^/3^A- 1 /m! 

m 

= e - A ''on e 7Pln(ro/r 1 ) ^ (/^A^ ) m /m! 

m 

p-'Wi „7Pln(ro/'*i) ( o/3 p A rori 



Using equation (A4), we have 

(^„n) = (ro/ri)"^, 

with 

$(p) = - 7 |p-(l-^)/(l- ) 9)]. 

Therefore 

€(p) 



For moments equation, we have 

tW 2 v \ 

= (r/r ) 



Therefore, 



and 



(5p 2p ) 

In — r . = K p In r + const, 

For scale-scale correlation, we have 

(5p p 5(f ) (W p W p ) 

ri ' r2 ~ WWK) ~ 

Using equation (A10), if keeping ri/r 2 to be constant, the rx-dependence of C p [ 
by 

C^=^2/ri)r«* ) - 2CW . 
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